Multiple-Source Multiple-Sink Maximum Flow in 
Directed Planar Graphs in 0(n L5 logn) Time 

Shay Mozes 
Brown University 

December 30, 2010 

Abstract 

We give an 0(n logn) algorithm that, given a directed planar graph 
with arc capacities, a set of source nodes and a set of sink nodes, finds a 
maximum flow from the sources to the sinks. 

1 Introduction 

In this paper we give an 0(n 15 log n)-time algorithm for the maximum flow 
problem in directed planar graphs with multiple sources and sinks (MSMS). 

• Input: 

— a directed planar embedded graph G with non-negative arc capacities 

— a set S of source nodes 

— a set T of sink nodes 

• Output: a feasible flow in G that maximizes the total flow into T. 

Several new results were recently obtained for this and related problems. 
A few months ago Borradaile and Wulff-Nilsen [1] and Klein and Mozes [5] 
independently presented two 0(n 15 logn)-time algorithms for the maximum 
flow problem in directed planar graphs with multiple sources and a single sink 
(MSSS). Very recently Nussbaum [7] presented a recursive 0(n 15 log 2 n)-time 
algorithm for MSMS. His algorithm uses, among other techniques, MSSS max 
flow computations. 

The proposed algorithm is similar to that of Nussbaum, but uses a simpler 
recursive approach. It uses a technique to redistribute excess flow among a cer- 
tain set of nodes called boundary nodes. This technique was used by Nussbaum 
for the case where the sources and sinks are separated by the boundary nodes. 
Here we prove that it is applicable in the general case. 
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2 preliminaries 



2.1 Jordan Separators for Embedded Planar Graphs 

Miller [6] gave a linear-time algorithm that, given a triangulated two-connected 
n-node planar embedded graph, finds a simple cycle separator consisting of at 
most 2\[2%Jn nodes, such that at most 2n/3 nodes are strictly enclosed by the 
cycle, and at most 2n/3 nodes are not enclosed. 

For an n-node planar embedded graph G that is not necessarily triangulated 
or two-connected, we define a Jordan separator to be a Jordan curve C that 
intersects the embedding of the graph only at nodes such that at most 2n/3 
nodes are strictly enclosed by the curve and at most 2n/3 nodes are not enclosed. 
The nodes intersected by the curve are called boundary nodes. To find a Jordan 
separator with at most 2^/2y/n boundary nodes, add artificial edges with zero 
capacity to triangulate the graph and make it two-connected without changing 
the maximum flow in the graph. Now apply Miller's algorithm. 

A cycle separator G separates the graph G into two subgraphs Gi , G2 called 
pieces. Gi is the embedded subgraph consisting of the nodes and edges enclosed 
by G, i.e. including the nodes intersected by C . Similarly, G2 is the subgraph 
consisting of the nodes and edges not strictly enclosed by G, i.e. again including 
the nodes intersected by C. 

2.2 Flow 

Let G be a directed graph with arc set A, node set V, source set S C V and 
sink set TCF-S. For notational simplicity, we assume here and henceforth 
that G has no parallel arcs and no self-loops. 

We associate with each arc a two darts d and c?', one in the direction of 
a and the other in the opposite direction. We say that those two darts are 
reverses of each other, and write d — rev(cf). Given an arc a with capacity c, 
the capacity associated with the drat d going in the direction of a is c. The 
capacity associated rev(d) is zero. 

A flow assignment /(•) is a real- valued function on darts that satisfies anti- 
symmetry: 

/(rev(rf)) = -f(d) 

A flow assignment /(•) is said to respect capacities if, for every dart d, f(d) < 
c(d). Such a flow assignment is also called a pseudoflow. 

For a given flow assignment /(•), the net inflow (or just inflow) of a node v 

is inflow/ (v) = Y,aeA-.hesid(a)=v f( a ) ~ EaeA:tail(a)=. /(»)• 

A flow assignment /(•) is said to obey conservation if for every node v (ji SUT, 
inflow/ (t>) = 0. A pseudoflow that obeys conservation is called a feasible 
flow. The value of a feasible flow /(•) is the sum of inflows at the sinks, 
X^tgT mn ow/ (t). A maximum flow is a feasible flow whose value is maximum. 

The residual graph of G with respect to a flow assignment /(•) is the graph G/ 
with the same arc-set, node-set, sources and sinks, and with capacity assignment 
c/(-) defined by Cf{d) = c(d) — f(d) for every dart d. 
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3 The Algorithm 



We present the algorithm as a recursive procedure with calls to the following 
subroutines: 

• MultiplcSourccSinglcSmkMaxFlow(G, S, t) - computes a maximum flow 
in G from source set S to sink node t. This can be implemented in 
0(n 15 logn) time using the algorithm in [5] or Q]. 

• SingleSourceMultiplcSinkMaxFlow(G, s, T) - computes a maximum flow 
in G from source node s to sink set T. This can be implemented using 
the multiple-source single-sink maximum flow subroutine. 

• SingleSourceSingleSinkLimitedMaxFlow(G, s, t, A) - computes a flow in 
G from source node s to sink node t, whose value is the minimum between 
A and the value of the maximum s-to-t flow. This can be implemented in 
linear time when s and t are incident to the same face [3J [3] . 

We omit discussion of the base case of the recursion (the case where the graph 
size is smaller than a certain constant). Each of the recursive calls operates on a 
subgraph of the original input graph. We assume one global flow assignment /(•) 
for the original input graph, and one global capacity assignment c(-). Whenever 
a subroutine is called, it takes as part of its input the current residual capacity 
function c/(-), computes a flow assignment /(•), and then updates the global 
flow assignment /(•) by f(d) := f(d) + f(d) for every dart in the subgraph. In 
the pseudocode, we do not explicitly mention /(•), c(-) er(-), c/(-), or The 
pseudocode for the algorithm is given below. 

The algorithm first finds a Jordan separator C with pieces G\ and G2. For 
each piece, it calls itself recursively (Line [3]), so that after the call there are no 
S-to-T residual paths within G^. It then pushes flow from the sources in Gj to 
the boundary nodes using a multiple-sources single-sink max flow computation 
(Line [5]), and similarly pushes flow from the boundary nodes to the sinks in G,; 
(Line After the first loop terminates, there are no S'-to-G residual paths, no 
G-to-T residual paths and no S-to-T residual paths in the entire graph. The 
resulting flow assignment is a pseudoflow rather than a feasible flow since it does 
not satisfy conservation at the boundary nodes. 

The algorithm then handles the boundary nodes one by one in cyclic or- 
der. Along the iterations of the second loop, we say that a boundary node pj 
is processed if the current value of the loop variable i is greater than j, and 
unprocessed otherwise. If the inflow at an unprocessed node pi is positive, the 
algorithm tries to resolve the violation of conservation at Pi by sending the ex- 
cess flow from pi to other unprocessed nodes. Similarly, if the inflow at pi is 
negative, the algorithm tries to send flow from other unprocessed nodes to pi. 
This approach for resolving excess flow in the boundary nodes was very recently 
used in [7] for the special case where all sources are in one piece and all sinks 
are in the other. We prove that the procedure works even in the general case 
where the sources and sinks are in both pieces. 
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Algorithm 1 MultipleSourccMultiplcSinkMaxFlow(graph G, sources S, sinks 



1: find a simple cycle separator G in G with pieces G\ and G2 
2: for i — 1, 2 do 

3: MultipleSourceMultipleSinkMaxFlow(Gi, SnG„ Tfl Gj) 

4: add to Gi artificial bi-directional arcs with infinite capacity between the 

boundary nodes and an artificial node v* embedded in the face of Gi 

resulting from the deletion of arcs of G not in Gi . 
5: MultipleSourceSmgleSinkMaxFlow(G l ,S' n G i: v*) 
6: SingleSourceMultipleSinkMaxFlow(Gi,?;*,T n G») 
7: remove v* and the artificial arcs from Gi 
8: Let the nodes of C be pi,p%, ■ ■ ■ ,Pk 
9: for i = 1, 2, . . . , k do 

10: Add infinite capacity artificial bi-directional arcs pi + ij>i + 2, pi+2Pi+3, • ■ •> 

Pk-lPk 

11: if inflow (pj) > then 

12: SingleSourceSingleSinkLimitedMaxFlow(G, pi : pi+\ , inflow (p^)) 
13: else 

14: SingleSourceSingleSinkLimitedMaxFlow(G,pi+i,Pi, |inflow(p,)|) 
15: remove artificial arcs 

16: push flow from boundary nodes with positive inflow back to sources and to 
boundary nodes with negative inflow from sinks 



After all boundary nodes are processed, the algorithm converts the pseud- 
oflow to a maximum feasible flow by sending any remaining excess flow from 
the boundary back to the sources and filling flow deficiencies by sending flow to 
the boundary from the sinks. This can be done in linear time by first canceling 
flow cycles using the technique of Kaplan and Nussbaum [J] , and then pushing 
the flow in topological sort order (cf. [5]) 

3.1 Correctness 

We will use the following two lemmas in the proof of correctness: 

Lemma 3.1 (suffix lemma) Let f be a flow with source set X. Let A, B be 
two disjoint sets of nodes. If there are no A-to-B residual paths and no X-to-B 
residual paths before f is pushed, then there are no A-to-B residual paths after 
f is pushed. 

Proof: / may be decomposed into a cyclic component (a circulation) and an 
acyclic component. Note that pushing a circulation does not change the amount 
of flow crossing any cut. This implies that if there were no A-to-B residual paths 
before / was pushed, then there are none after just the cyclic component of / 
is pushed. It therefore suffices to show the lemma for an acyclic flow /. 
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Suppose, for the sake of contradiction, that there exists a residual a-to-6 
path P after / is pushed for some a £ A and b £ B. Let P' be the maximal 
suffix of P that was residual before the push. That is, the arc e of P whose 
head is start(P') was non-residual before the push, and /(rev(e)) > 0. The 
fact that /(rev(e)) > implies that before / was pushed there was a residual 
path Q from some node x £ X to head(e). Therefore, the concatenation of Q 
and P' was a residual ir-to-& residual path before the push, a contradiction. QED 

Lemma 3.2 (•prefix lemma) Let f be a flow with sink set X. Let A, B be two 
disjoint sets of nodes. If there are no A-to-B residual paths and no A-to-X 
residual paths before f is pushed, then there are no A-to-B residual paths after 
f is pushed. 

Proof: Similarly to the proof of Lemma 13.11 it suffices to prove the lemma 
for an acyclic flow /. Suppose, for the sake of contradiction, that there exists a 
residual a-to-6 path P after / is pushed for some a £ A and b £ B. Let P' be 
the maximal prefix of P that was residual before the push. That is, the arc e 
of P whose tail is end(P') was non-residual before the push, and /(rev(e)) > 0. 
The fact that /(rev(e)) > implies that before / was pushed there was a resid- 
ual path Q from tail(e) to some node x £ X . Therefore, the concatenation of 
P' and Q was a residual a-to-X residual path before the push, a contradiction. 
QED 

As Nussbaum points out [7] , we may assume, without loss of generality, that 
no sources or sinks belong to C. Otherwise we may replace each such terminal 
(i.e., either a source or a sink) v with a new terminal v' embedded in a face to 
which v is adjacent, connect v' to v and designate v' as the terminal instead of 
v. 

Lemma 3.3 After Line\3\is executed for piece Gi, there is no S-to-T residual 
path in Gi. 

Proof: By maximality of the flow pushed in Line [3] QED 

Lemma 3.4 After Line\^ is executed for piece Gi, there is no: 

1. S-to-T residual path in Gi 

2. S-to-C residual path in Gi 

Proof: Item [2] follows from the maximality of the flow pushed in Line[5j ItcmQ] 
follows by applying the suffix lemma with A = SC\G l , B = TC\Gi, X = SnG^ 
and / the flow pushed in Line [5] QED 
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Lemma 3.5 After Line\d\is executed for piece Gi, there is no: 

1. S-to-T residual path in Gi 

2. S-to-C residual path in Gi 

3. C -to-T residual path in Gi 

Proof: Item [3] is immediate from the maximality of the flow pushed in Line[6l 
Item [T] follows by applying the prefix lemma with A = S n Gi, B = T n Gi, 
X = T n Gi, and / the flow pushed in Line [fO Item [2] follows by applying the 
prefix lemma with A = S H Gi , B = C, X = T n Gi , and / the flow pushed in 
LineH QED 

An immediate corollary of Lemma 13.51 is 
Corollary 3.6 Immediately after the first loop terminates there is no: 

1. S-to-T residual path in G 

2. S-to-C residual path in G 

3. C-to-T residual path in G 

Lemma 3.7 The following invariants are preserved throughout the execution of 
the second loop 

1. There is no S-to-T residual paths in G. 

2. There is no residual S-to-C path nor residual C-to-T path in G. 

3. If a processed node pj has positive inflow, there is no residual path from 
Pj to the as-yet-unprocessed nodes. If Pj has negative inflow, there is no 
residual path to it from the as-yet-unprocessed nodes. 

4- There is no residual path from a processed node with positive inflow to a 
processed node with negative inflow. 

Proof: By induction on the number of iterations i of the loop (i.e., the num- 
ber of processed nodes). By corollary 13.61 the first two invariants are satisfied 
immediately before the second loop is executed. The last two invariants are 
trivially satisfied since at that time there are no processed nodes. 

Assume the invariants hold up until the beginning of the i th iteration. Sup- 
pose that pi has positive inflow at the beginning of the iteration (the case of 
negative inflow is similar). At the end of the iteration pi is processed and we 
need to show that the invariants still hold. 

1. Invariant [T] holds by invoking the suffix lemma with A = S, B = T, 
X = {pi}, and / the flow pushed from pi in LincfT^l 
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2. There are no residual S-to-C paths by invoking the prefix lemma with 
A = S, B = C, X = {pj : j > i} and / the flow pushed from pi in 
Line [HJ There are no residual C-to-T paths by invoking the suffix lemma 
with A = C, B — T, X — {pi} and / the flow pushed from pi in Line[T2l 

3. Since pi had positive inflow at the beginning of the iteration, and the flow 
pushed in Line Q2] is limited, if pi has non-zero inflow at the end of the 
iteration it must be positive, and the flow pushed was in fact a maximum 
flow from pi to {pj : j > i}. Invariant [3] holds for pi by maximality of the 
flow pushed in Line 1121 The invariant holds for processed nodes pj with 
j < i by invoking the prefix lemma with A = {pj : j < i}, B — {pj : j > i}, 
X = {pj : j > i}, and / the flow pushed from pi in Line 1121 

4. Invariant 2] holds for {pj : j < i} by invoking the prefix lemma with 
A = {pj : j < i, inflow(pj) > 0}, B = {pj : j < i, inflow (pj) < 0}, 
X = {pj : j > i}, and / the flow pushed from pi in Line ll2l The invariant 
holds for pi by invoking the suffix lemma with A = {pA-, B = {pj : j < 
i, inflow(pj) < 0}, X = {pi} and / the flow pushed from pi in Line 1121 

QED 

Theorem 3.8 The flow computed by the algorithm is a maximum feasible flow. 

Proof: The flow pushed by the algorithm originates only at sources and bound- 
ary nodes and terminates only at sinks and boundary nodes. Therefore, sources, 
sinks and boundary nodes are the only nodes whose inflow might be non-zero. 
Since line [TH] makes the inflow at all boundary nodes zero, the flow assignment 
upon termination is a feasible flow. It remains to show that upon termination 
there is no residual S'-to-T path. Let C+ (C_) be the set of nodes with posi- 
tive (negative) inflow just before Line [TB] is executed. Let / + (/_) be the flow 
pushed back from C + to S (from T to C_ ) in Line [T5] Consider first pushing 
back / + . By Lemma 13. 71 we may invoke the suffix lemma with A — S U C_ , 
B = T, X = C + and / = / + to show there are no S'-to-T residual paths nor 
C_-to-T residual paths after / + is pushed. Similarly, invoking the suffix lemma 
with A = S, B = C-, X — C+ and / = /+ shows there is no S-to-C- residual 
path after /+ is pushed. Next, consider pushing /_. Invoking the prefix lemma 
with A — S, B = T, X = C- and / = /_ shows there are no residual S-to-T 
paths after Line [1^] is executed. QED 



3.2 Running Time 

The algorithm performs one recursive call per piece. In addition, it performs 
two multiple-source single-sink max flow computations per piece, which take 
0(|Gi| logn) time. Processing the boundary nodes takes 0(n) per node since 
the maximum flow computed is between a source and a sink on the same face. 
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Therefore, the non-recursive part takes 0(n 15 logra). Since the size of the two 
pieces (including boundary nodes) are at most On + and (1 — 6)n + 

for some 1/3 < 9 < 2/3, the total running time is thus bounded by 

T(n) < cin 15 logn+ max \t(9u + VSn) + T((l - 6)n + VSn)} (1) 

l/3<$<2/3 v J 

for some constant c\. 

Lemma 3.9 T(n) = 0(n 15 logn) 



Proof: We prove, by induction on n, that there exists a constant c such that 
for sufficiently large n, T(ri) < en 15 log n. We choose c sufficiently large so that 
the base of the induction, where n is some constant, holds. For the inductive 
step, using the inductive hypothesis in Eq. [T]we get 



T(n) < n L b \ogn 



ci 



max 

l/3<#<2/3 



/ 8/n 



1.5 



By convexity, the maximum is attained at the extreme values. 
Since (1/3) 15 + (2/3) 15 = 0.7376 . . . , for sufficiently large n, (o + 



1.5 



/ \ 1-5 

1 1 — 9 + y 8/n J < 0.74. Therefore, c can be chosen sufficiently large so that 



C\ + 0.7 Ac < c, proving the lemma. QED 



We note that the running-time bottleneck are the multiple-source single-sink 
max flow computations. An 0(n . ) bound on MSSS would imply, by the above 
proof an 0(n 15 ) bound for multiple-sources multiple-sinks maximum flow as 
well. 
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